Transition times in the Landau-Zener model 
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This paper presents analytic formulas for various transition times in the Landau-Zener model. 
Considerable differences are found between the transition times in the diabatic and adiabatic bases, 
and between the jump time (the time for which the transition probability rises to the region of its 
asymptotic value) and the relaxation time (the characteristic damping time of the oscillations which 
appear in the transition probability after the crossing) . These transition times have been calculated 
by using the exact values of the transition probabilities and their derivatives at the crossing point and 
approximations to the time evolutions of the transition probabilities in the diabatic basis, derived 
earlier [N. V. Vitanov and B. M. Garraway, Phys. Rev. A 53, 4288 (1996)], and similar results in 
the adiabatic basis, derived in the present paper. 
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I. INTRODUCTION 



The Landau-Zener (LZ) model JJg] has long ago be- 
come a standard notion in quantum physics. It provides 
the probability of transition between two quantum states 
coupled by an external field of a constant amplitude and 
a time-dependent frequency which passes through reso- 
nance with the transition frequency. This level cross- 
ing, seen in the diabatic basis (the basis of the two bare 
states - the eigenstates of the Hamiltonian in the ab- 
sence of interaction), appears as an avoided crossing in 
the adiabatic basis (the basis comprising the two eigen- 
states of the Hamiltonian in the presence of interaction) . 
Cases of level crossings and avoided crossings can be met 
in a number of areas in physics, such as quantum optics, 
magnetic resonance, atomic collisions, solid-state physics, 
atom-surface scattering, molecular physics, and nuclear 
physics. The Landau-Zener model is a reliable qualita- 
tive (and often even quantative) tool for describing and 
understanding such phenomena. 

Along with the transition probability, it is often nec- 
essary to know the time for which the transition occurs 
- the transition time. For example, in the case of multi- 
ple level crossings (or avoided crossings) it is essential to 
know if the transition is completed before the next cross- 
ing. Moreover, even in the case of a single crossing the 
transition time is an important parameter because the 
actual coupling and detuning never match the LZ ones 
exactly. For instance, the LZ formula applies when the 
transitions take place in a narrow time interval around 
the crossing, provided the actual detuning changes nearly 
linearly with time and the coupling is nearly constant in 
the vicinity of the crossing. However, when transitions 
can occur far from the crossing, where the actual coupling 
usually is not constant and the detuning is not linear in 
time, the usage of the LZ formula may be incorrect. 

It is far from obvious what is the transition time in 
the diabatic basis because there the coupling Q is con- 
stant and lasts from — oo to +oo. Since the detuning, 
being a linear function of time, A = f3 2 t, diverges at ±oo, 



the transition probability is well defined, but there is no 
apparent time scale at which the transition takes place. 
It looks considerably easier to determine the transition 
time r a in the adiabatic basis because there the (non- 
adiabatic) coupling is a Lorentzian function of time with 
a width of 20,/ 1 and the energies of the two adiabatic 
states have an avoided crossing with the same duration; 
hence, the transition time is expected to be r Q w 2f2//3 2 . 
This deceptively obvious conclusion turns out to be only 
partially correct. The problem here arises from the fact 
that the non-adiabatic coupling vanishes too slowly and 
also, the eigenenergy gap diverges too slowly. 

The scaling properties of the LZ transition time t c i in 
the diabatic basis have been studied by Mullen et al [[|, 
who have found that for large f2, is proportional to 
Q,, while for small f2, is constant. These authors have 
used two different approaches - an "internal clock" , based 
on perturbative calculation of the transition probability 
time evolution, and an "external clock" , based on identi- 
fying a characteristic frequency appearing in the response 
of the two-state system to a harmonic perturbation. 

In the present paper, I derive analytical estimates for 
the LZ transition times by using some exact and approxi- 
mate results for the transition probability in the diabatic 
basis, derived in and similar results in the adiabatic 
basis derived here. Thus, this paper provides not only 
the scaling properties but also explicit formulas for the 
transition times in both the diabatic and adiabatic bases. 
In view of the numerous exact and approximate results 
for the transition probabilities, it turns out much easier 
to calculate the transition time than to define it. I distin- 
guish two kinds of transition times. The jump time r Jump 
is the time for which the transition probability rises to 
the region of its asymptotic value P(oo) (the exact def- 
inition is given below). The relaxation time r rclax is the 
time for which the amplitude of the oscillations, which 
appear in the transition probability after the crossing, is 
damped to a sufficiently small value, eP(oo), (e <C 1). 

The paper is organized as follows. The basic equations 
and definitions are given in Sec. ||. The transition times 



1 



[V. The conclusions are 



in the diabatic basis are calculated in Sec. Ill and those 
in the adiabatic basis in Sec. 
summarized in Sec. |v|. 



II. BASIC EQUATIONS AND DEFINITIONS 

A. Diabatic basis 

The time evolution of a coherently driven two-state 
quantum system is described by the Schrodingcr equation 
(in units h = 1) J§ 



d 


"<*(*)" 




i — 

dt 


_c 2 (t)_ 





-A(t) Q(t) 
Q(t) A(t) 



Cl (t) 

C2(t) 



(1) 



where c\(t) and c 2 {t) are the probability amplitudes of 
states tpi and tp 2 , &(t) is the coupling (assumed real) 
between the two states and A(t) is a half of the difference 
between the system transition frequency and the field 
frequency. In the Landau-Zener model, we have 



Q(i) = const, A(t) 



(2) 



where the coupling £!(i) is supposed to last from t — > — oo 
to t — > +oo. Without loss of generality the slope 1 of 
the detuning as well as the real constants f2 and j3 are 
assumed positive. Both £1 and (3 have the dimension of 
frequency. Following the notation of |4| , I introduce the 
scaled dimcnsionless time r and coupling u>, 
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Provided the system has been initially in state "01: the 
probability of transition to state %p 2 at time r, P&{t) = 
|c 2 (r)| 2 ,is||-§ 

Pd(r) = ye-» 2 / 4 \D_ w/2 (rV2 t 
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(4) 



where D v {z) is the parabolic cylinder function QJ^]. Here 
the subscript "<f ' indicates the diabatic basis. 



B. Adiabatic basis 



The two adiabatic states (defined as the instanta- 
neous eigenstates of the Hamiltonian) are given by 
y>i(r) = f/'i cos ??(t) — ip 2 sin $ (r), <p 2 {r) — V>isin$(r) + 
-02 cos i9(r), where the angle $(t) is defined as 



tan2??(r) 



Q(r) 
A(r) 



r 



(0 ^ 0(t) ^ tt/2). (5) 



The eigenvalues of states ip% and ^2 are given by — Q.q 
and +Oo, respectively, and the (non-adiabatic) coupling 
between them by i?'(t) = di9(r)/dr, where 



J2 (i) = v^Ti 



tf'(r) 



2(r 2 +w 2 ) 



(6) 



The condition for adiabatic evolution is \$'\ <C £lo, which 
requires that uj 2 ^> 1. Hence, the coupling uj = Q/f3 plays 
also the role of the adiabaticity parameter. 

The transition probability Pb.{t) between the two adi- 
abatic states can be obtained from the relation between 
the evolution matrices in the adiabatic and diabatic bases 
by using the exact diabatic evolution matrix given in 0] 
and some properties of D v (z) M. The result is 
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/2 (rV^e^/ 4 ) cosz9(r) 
/4 D_ 1+tLu2/2 (rV2e 3 ^ 4 ) sintf(r) 



(7) 



In the figures below, both the diabatic and adiabatic 
transition probabilities are calculated by the highly accu- 
rate numerical method desribed in the Appendix, rather 
than from Eqs. (Q) and (J?]). These equations are used for 
analytic derivation of the transition times. 



C. Some exact values of the transition probabilities 

For the calculation of the transition times, we need a 
few values of P^t) and -P a (r), easily obtained from Eqs. 
([|) and (0). By taking the limit r — ► oo, we recover the 
well-known LZ probabilities 



Pd(oo) = 1 



^a(oo) 



(8) 
(9) 



By using the power series |4|]8|,|9| of D v (z) we can expand 
Pd(r) and P&(t) in terms of r, which enables us to find 
the values of Pi(r) and P a (r) and their derivatives at 
t = 0. We need the following values 



Pd(0) 
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P^(0) = oj\/l - e-™ 2 cos\, 
P^'(0) = 2uj 2 e-™ 2/2 , 



and 



fa(0) 



i-VT- 



sinx 



P"(0) = VI - e"™ 2 



smx 
2lo 2 



cosx 



where the angle x( w ) is defined by 

= J + argr (i - ii W 2 ) - argr (l 



(10a) 

(10b) 
(10c) 

(11a) 
(lib) 
(11c) 

. (12) 
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It is a monotonically increasing function of uj. For small 
or large uj, x{uj) behaves as 



\<-=^ + ^ 2 -#" 6 +^ 10 ). (- 2 «1), (13a) 
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x(-)=2-2^-3^ + ^ U )»(^» 1 ). ( 13b ) 
where £(z) is the Riemann's zeta function |H|. 



III. TRANSITION TIMES IN THE DIABATIC 
BASIS 

A. Time evolution of the transition probability 

The time evolution of the diabatic transition probabil- 
ity Pd (t) is shown in Fig. [l] for six values of the coupling 
uj - from 0.03 (small) to 10 (large). The evolution shows 
two characteristic time regions. The first one is around 
the crossing (r = 0), where Pa(t) rises from zero to about 
its asymptotic value Pd(oo); this region determines the 
jump time r^ ump . This jump is followed by a region, 
where Pa{t) oscillates around the value Pd(oo) (for large 
uj, these oscillations become invisible); this region deter- 
mines the relaxation time r™ lax . 

The time evolution of Pd(i") before and after the cross- 
ing is well approximated by the formulas B 



P d (r<0)« | + 
P d (r>0)« |+ (§- 



2Vr 2 + uj 2 ' 



Vt 2 + UJ 2 



\/t 2 + uj- 



(14a) 



: cos£(t), (14b) 



valid for r 2 + uj 2 > 1. The phase £(t) is given by 

,2 

— L i.j 2 In 

s/2 



e(r)=-^+c 2 ln 



4= ( r + Vt 2 + uj' 



+ry/r 2 + uj 2 + J + argT (l - \iuj 2 ) . (15) 



B. Jump time 

The attempt to define the jump time in a simple and 
unambiguous manner quickly comes across some difficul- 
ties. First of all, it is not obvious how to define the initial 
time of the transition, because Pa(t) is nonzero for any 
finite time. A reasonable choice seems to be the time 
when the rising transition probability P<i(t) first equals 
eP d (oo), where e is a suitably chosen small positive num- 
ber (e <C 1). It is even less obvious how to define the final 
time of the transition in a way that applies to both small 
and large coupling. One possible choice, used by Lim and 
Berry fefl , is to take the time at which the upper envelope 
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FIG. 1. The time evolution of the diabatic transition prob- 
ability -Pd(i~) (thick curve) for six values of the coupling uj - 
from uj = 0.03 (small) to u — 10 (large). In each figure, the 
horizontal dashed line shows the asymptotic value Pd(oo), Eq. 
(^) [for uj — 3 and 10, it coincides with the axis -Pd(i") = 1]. 
The two vertical dashed lines display the jump region, whose 
beginning is defined by the crossing of the tangent to P&{t) 
at t = (shown by a slanted solid line) with the zero line, 
whereas its end is defined by the crossing of the tangent with 
the Pd(oo) line. 



of t he os cillations [obtained by setting cos £(t) = — 1 in 
Eq. ( |l4b| )] touches unity. However, this idea does not ap- 
ply for small uj, because then the upper envelope never 
touches unity. Another possibility is to take the time at 
which Pi(t) crosses its asymptotic value Pd(oo) for the 
first time. However, this is appropriate for small uj only, 
because for large uj (when the oscillations are strongly 
damped) , the crossover takes place at exponentially large 
times, although Pd(7") comes very close to Pd(oo) much 
earlier. Alternatively, one can take the time when Pd(7") 
first equals the value (1 — e)Pj(oo), which is reasonable 
and consistent with the e-definition of the initial time 
discussed above, but leads to a rather complicated ex- 
pression. I propose here a more elegant and simple solu- 
tion, which applies to any uj and provides similar results 
as the e-approach. I define the jump time as 



jump _ 



Pd(oo) 
^(0) ' 



(16) 



This definition is based upon the geometrical meaning of 
the derivative as the slope of the function at the point of 
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FIG. 2. The jump time (^) and the relaxation time ( |is| ) 
(for e — 0.1) of the diabatic transition probability, plotted 
against the dimensionless coupling ui = £1/13. 



calculation. It provides good results when applied to the 
most frequently used smooth functions that rise mono- 
tonically from to 1, e.g. f(x) — |(1 + tanhx) and 
f(x) = i(l + xj\l ' x 2 + 1), typically providing the inter- 
val where the function rises from about 0.10 — 0.15 to 
about 0.85 — 0.90. By using the e xact values of Pa(oo) 



and P^(0) from Eqs. (Jq) and (tOq), we obtain 



jump _ 



uj cosx(oj) 



It follows from Eqs. (13a) and (13b) that 



lump 
jump 



V^F, {lo 2 « 1), 
2ui, {uj 2 > 1). 



(17) 



(18a) 
(18b) 



In other words, the jump time is proportional to uj at 
large uj while it is nearly constant for small uj. Thus, Eqs. 

P, and (18b) confirm the scaling properties found in 
>r the extreme cases of small and large uj. 
This behavior of the jump ti me f or sm all a nd large 
ui can be explained from Eqs. ( |l4a| ) and (14b), which 
provide Pd(r < 0) and Pd{r > 0), and from the Taylor 
expansion of Pd(j~) around r = 0, obtained by using the 
derivatives ([l0|). It can easily be shown that for large 
w, P&{t) depends on the ratio t/uj only which means 
that r d ump oc uj. For small uj, the normalized transition 
probability Pd(r)/Pd(oo) depends on r only, which can 
indeed be seen in Fig. | for uj = 0.03, 0.1 and 0.3; hence, 
r d ump must not depend on uj. 
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FIG. 3. The adiabatic transition probability evolution 
Ps.(t) (solid curve) plotted against the scaled time t/uj for 
the same six values of the coupling uj as in Fig. |l| - from 
uj = 0.03 (small) to u = 10 (large). The two vertical dashed 
lines in each figure display the interval — u) ^ r ^ uj, which 
provides approximately the jump region for uj < 1. 



I define the relaxation time r, 



relax 



as the time it takes 



to damp the oscillation amplitude t o the (small) value 
ePd(oo), where £<?;!. By using Eq. ( [l4b] ), we find 



(19) 



-1 



The square root is real only for ttuj 2 ^ ln(l/e 2 + 1). This 
inequality imposes an upper limit of uj, above which the 
oscillation amplitude is never larger than ePd(oo). For 
e = 0.1, this limit is uj < 1.21. 

The diabatic jump time ( p7| ) and the relaxation time 
(|l9| ) (with e = 0.1) are displayed in Fig. |[ 



IV. TRANSITION TIMES IN THE ADIABATIC 
BASIS 

A. Time evolution of the transition probability 



C. Relaxation time 



As Eq. ( p_4b| ) shows, the amplitude of the oscillations 
in Pd(r > 0) vanishes as r _1 at large positive times. 



The time evolution of the adiabatic transition proba- 
bility Ps,(t) is shown in Fig. || for the same six values 
of the coupling uj as in Fig. [l| - from 0.03 (small) to 10 
(large). There are two distinct types of evolution. For 
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small ui, -P a (r) behaves as the diabatic transition prob- 
ability Pd( T ) m Fig. [A For large lo, P a (r) rises from 
zero at — oo to its maximum near the crossing (r = 0) 
and then decreases to its exponentially small asymptotic 
value P a (oo) = e _7rw (Hi. The small-u) case can be 
treated in the same manner as in the diabatic basis, while 
the large-w case requires a more carefull analysis. 

The time evolution of the adiabatic transition proba- 
bility before and after the crossing is approximated by 
the formulas 



P a (r < 0) 



16(r 2 + lo 2 ) 3 
P a (r > 0) sa e~™ 2 + ( 1 - 2. 



(20a) 



16 (t 2 +lo 2 Y 



2(r 2 +w 2 )2 



■sin£(r), (20b) 



valid for r 2 + lo 2 > 1. The phase £(r) is given by Eq. 
@. Equations @, which are new, can be derived from 
Eq. (0) in the same manner as Eqs. jl4| ) have been de- 
rived from Eq. (^) in 0), but by keeping more terms 
in the large- argument- and-large- order asymptotic expan- 
sions |ll]] of the parabolic cylinder functions in Eq. (R) . 



B. Jump time 

1. Small uj 

For small lo, the transition probability evolution re- 
sembles that in the diabatic basis. Hence, I define the 
jump time r| ump in the same way as Eq. (fL6h, 



r jump 



■P.(oo) 



By using Eqs. and ( lib ) we find that 
Ti ump = 2uje~™ 2/2 ~ °- • ' ' 2 



2w, 



(^ < 1). 



(21) 



(22) 



Hence, for small w, the jump time in the adiabatic basis 
is proportional to w, as expected. 



2. Large u> 

For large lo, I define the initial time of the transition 
as the time -rh™ 15 ' 1 < at which P a (r) = eP a (oo), where 
e is a suitably chosen small number. It follows from Eq. 



(20a) that 
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(23) 
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FIG. 4. The time evolution of the adiabatic transition 
probability P a (r) for uj — 2. The upper figure displays the 
region around the crossing (r = 0). The lower figure gives 
an expanded view of the region after the crossing, which 
shows the upper and lower envelopes of the oscillations and 
the non-oscillatory part of Pa(i") (short-line dashed curves). 
The horizontal solid line depicts the asymptotic value P a (oo), 
whereas the two horizontal dashed lines above and below it 
show the values (1 + e)P a (oo) and (1 — e)P a (oo), respectively, 
with e — 0.1. The two vertical dashed lines show the fi- 
nal time of the jump r a ump ' f [defined by the crossing of the 
non-oscillatory part of Pa.{r) with the (1 + e)P a (oo)-line] and 
the relaxation time r a elax [defined by the time when the os- 
cillation amplitude gets smaller than eP a (oo)]. 

To define the final time of the jump r^ ump,f > 0, I 
first remark that as r increases after the crossing , the 
non-oscillatory part of the transition probability fl20b| ) 
approaches the asymptotic value P a (oo) from above [be- 
cause f or > 2, i. e., for lo > 0.47, the second term 
in Eq. (20b) is positive]. Hence, I define rj ump,f as the 



time at which the non-oscillatory part of Pa.{j) is equal to 
(1 + e)P a (oo). An illustration of this definition is shown 
in Fig. ^. A simple calculation gives 



T jump,f _ UJA 



1 

2\ 3 



The total jump time is 



T jump _ r jump,f _ T jump,i 



(24) 



(25) 



For (J 2 3> 1, we have 
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2 /6 



(lo 2 » 1). 



(26) 



Thus, for large lo, the transition time in the adiabatic 
basis is not proportional to lo, but it rather increases 
exponentially. This behavior can be explained by the fact 
that for large lo, P a (oo) (= e _7r " ) is exponentially small 
and then the population changes in the slowly vanishing 
wings of the non-adiabatic coupling i?'(r) [see Eq. (j^)] 
are non- negligible compared to P a (oo). 



C. Relaxation time 

As Eq. ( p0b| ) shows, at large positive times the ampli- 
tude of the oscillations, that appear in P a (r) after the 
crossing, vanishes as t~ 3 . The relaxation time r^ olax is 
defined in the same way as in the diabatic basis - as the 
time it takes to damp the oscillatio n am plitude to the 
(small) value eP a (oo). By using Eq. fl20b|) , we find 



rf lax « w 



\ 



1 



4e 2 w 4 



For small and large lo, this equation reduces to 



.rcla 



.rcla 



(lo 2 » 1). 



(27) 

(28a) 
(28b) 



A comparison of Eqs. ( ft9| ) and ( [28a| ) shows that for 
small lo, Tj olax 3> r™ lax . This is explained by the fact 
that the oscillation amplitude of P a (r) vanishes as t~ 3 , 
i. e., much faster than that of Pd(T) which vanishes as 
t _1 . In contrast, for large lo, we have r™ lax 3> T^ clax « 
0, which follows from the fact that the reference value 
in the diabatic basis is Pd(oo) = 1 — e~ nuJ w 1, while 
the reference value in the adiabatic basis is P a (oo) = 

The adiabatic jump time T^ mp and the relaxation time 
T a' clax (for e = 0.1) are display ed in Fig. ||. Note that, 
as follows from Eqs. ( l2u| ) and ( 28b ), the ratio between 
the jump and relaxation times at large lo is constant, 

T jump^ T reiax ^ ^6 g )e, i. e., they are almost equal for 
e = 0.1. In the evolution picture (Fig. ||), however, the 
relaxation ends later than the jump, because the jump 
time is calculated from r = r£ lmp l < 0, while the relax- 
ation time is calculated from the crossing (r = 0). 



V. SUMMARY OF THE RESULTS AND 
CONCLUSIONS 

In the present paper, I have calculated various tran- 
sition times in the Landau-Zener model. I have empha- 
sized the differences between the transition times in the 
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FIG. 5. The jump and relaxation times of the adiabatic 
transition probability, plotted against the dimensionless cou- 
pling lo = S7//3. The small-u; jump time (^) is shown by the 
solid line in the range ^ lo ^ 1, whereas the large-oj jump 
time (^B|) is shown by the solid curve for lo ^ 0.5. Note that 
the small-aj and large-tj formulas agree well for 0.5 < lo < 1. 
The relaxation time ( p7[ ) is plotted by a dashed curve. 



diabatic and adiabatic bases, and between the jump time 
and the relaxation time. The jump time is the time for 
which the transition probability rises to the region of its 
asymptotic value P(oo). The relaxation time is the time 
for which the amplitude of the oscillations, which appear 
in the transition probability after the crossing, is damped 
below the (small) value eP(oo), (e <C 1). These transi- 
tion times have been calculated by using the exact values 
of the transition probabilities and their derivatives at the 
crossing point as well as approximations to the transition 
probabilities evolutions derived in [Q and here. 

The results for the jump time in the diabatic basis 
j p con g rm the scaling properties found by Mullen et 
al p| in the limits of small and large coupling lo, i. e., 



for large lo, -H ump is proportional to lo, whereas for small 
lo, <rJ ump is constant. The jump time in the adiabatic 
basis T-j ump has a rather different dependence on lo. The 
seemingly obvious conclusion that rJ ump should be pro- 
portional to lo (because the non-adiabatic coupling is a 
Lorentzian function of time with a width of 2lo and the 
energies of the two adiabatic states have an avoided cross- 
ing with the same duration) turns out to be correct for 
small lo only, while for large lo, r^ ump increases exponen- 
tially. 

The relaxation times in the two bases, r| and r™ lax , 
also show rather different dependences on lo. The dia- 
batic relaxation time r™ lax is a decreasing function of lo 
and it vanishes above certain lo (« 1.2), whereas the adia- 
batic relaxation time r™ lax is an exponentially increasing 
function of lo. 

It should be pointed out that the transition times ob- 
tained in this work refer to the transition probabilities 
Pd(r) and P a (t). These may differ from the transition 
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times for the probabilities of no transition, 1 — P<i{t) and 
1 — P a (r), particularly those times which are linked to 
the values of the probabilities at r — > oo. 

The present paper has dealt with the transition times 
in the diabatic and adiabatic bases only, which are the 
most frequently used bases in practical applications of the 
LZ model. It has been shown by Berry and Lim fj|Jl2[ 
in their superadiabatic treatment of quantum evolution 
that the transition time is shortest and the oscillations 
in the corresponding transition probability are minimal 
in the optimal superadiabatic basis. 

In conclusion, the transition times obtained in this pa- 
per provide simple criteria for estimating the applicabil- 
ity of the Landau-Zener model to various cases of level 
crossings and avoided crossings. 



APPENDIX A: NUMERICAL INTEGRATION OF 
THE LANDAU-ZENER PROBLEM 

1. Diabatic basis 

Since in the LZ model (|^) the coupling does not van- 
ish at infinity and the detuning approaches infinity very 
slowly, the numerical integration of Eq. (|l|) is not a triv- 
ial problem, particularly when high accuracy is required. 
The straightforward way of integrating Eq. ([!]) is to start 
at a certain large negative time and propagate the so- 
lution towards +oo. However, a finite initial time t; 

generates spurious oscillations with an amplitude pro- 

_i 

portional to (r 2 + uj 2 ) 2 0] and one has to take a very 
large Tj in order to achieve a good accuracy in PdC?"), 
which is very expensive in terms of computation time. 
An alternative and much more efficient solution to this 
problem has been proposed in Q, which is summarized 
here for the reader's convenience. The transition prob- 
ability is derived from the equation for the population 
inversion w^(t) = 2P d {r) — 1 (derived from the optical 
Bloch equations p)), 



in 

TW d 



ll 

w d 



-At{u 2 + T 2 )w' d - Auj 2 w d = 0, (Al) 

rather than from Eq. ([!]) . The integration starts at r = 
and the solution is propagated towards the desired (pos- 
itive or negative) time. The initial conditions are found 
by identifying the terms in the Taylor expansion of Pd(i~) 
around r = [obtained by using the power series expan- 
sions of the parabolic cylinder functions in Eq. (^)] with 
the derivatives of Pd{T~) at r = 0. The initial values 
needed to start a Runge-Kutta algorithm, are 



2 /2 



w d (0) = -e 
w' d (0) = 2uj\J\ - e-™ 2 cosx 
w' d \0) = 4w 2 e-™ 2 / 2 



(0) = 4a; V 1 - e-™ 2 (sinx - 2oj 2 cosx) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



2. Adiabatic basis 

A similar numerical method, which is new and comple- 
ments the one described above for the diabatic basis [Q, 
can be used to obtain the transition probability P&(t) in 
the adiabatic basis and it has similar advantages. It turns 
out convenient to use the angle $ = \ arctan(o;/T), rather 
than the time r, as an independent variable. The equa- 
tion for the population inversion Wa[$(r)] = w a (r) = 
2-P a (i") — 1 has the form 



W'" + 6cot2tf W" 



4[4w 4 (cot 2 2t9 + l) 3 + 

+24cot2tf W a = 0, (A3) 



where a prime now means d/ctd. The initial values of 
W a (t?) and its derivatives at # = 7r/4 (i. e., at r = 0), 
needed to start a Runge-Kutta algorithm, are 



W a (7r/4) = -Vl - e" 7 ™ 2 sinx, (A4a) 

I<(7r/4) = -2e-™ 2 / 2 , (A4b) 

H^'(tt/4) = Wl-e-™ 2 (sinx - 2w 2 cosx) , (A4c) 

W a "(7r/4) = 8 (4w 4 + 1) e~™'/ 2 . (A4d) 

with x( w ) given by Eq. jl2}). 
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